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Abstract 

We discuss a recently demonstrated type of microwave spectroscopy of trapped 
ultra-cold atoms known as "echo spectroscopy" [M.F. Andersen et. al., Phys. Rev. 
Lett., in press (2002)]. Echo spectroscopy can serve as an extremely sensitive ex- 
perimental tool for investigating quantum dynamics of trapped atoms even when 
a large number of states are thermally populated. We show numerical results for 
the stability of eigenstates of an atom-optics billiard of the Bunimovich type, and 
discuss its behavior under different types of perturbations. Finally, we propose to 
use special geometrical constructions to make a dephasing free dipole trap. 
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1 Introduction 

The subject of decoherence of a quantum system has witnessed a renewed 
interest in the last years. In addition to being central to our understanding 
of the transition from quantum to classical physics[l], it is of outmost impor- 
tance in the blooming field of Quantum Information [2], where decoherence or 
dephasing represent loss of information. 

A quantum superposition state "collapses" due to the dissipative interaction 
with the environment (e.g. in the case of trapped atoms, the coupling to the 
electro-magnetic vacuum leads to spontaneous scattering of photons). This 
interaction is irreversible in nature. On the other hand, the macroscopic (en- 
semble averaged) response of a quantum system prepared in a superposition 
state decays due to the dephasing resulting from local variations in the evo- 
lution of the system (e.g. inhomogeneous broadening). By stimulating an ef- 
fective "time reversal", revivals of coherence in dephased systems (at least 
partial ones) have been reported for spin echoes [3] and photon echoes [4,5], 
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and more recently for a motional wave packet echo using ultra cold atoms in 
a one-dimensional optical lattice [6] . 

In this paper we discuss the technique of "echo spectroscopy" on ultra cold 
atoms trapped in an optical dipole trap, which utilizes a micro wave (MW) 
analog to the photon echo [7]. We show that this technique can be used to study 
the quantum dynamics of trapped atoms, even when many states arc thermally 
populated. We show that the stability of eigenstates under perturbations of 
the Hamiltonian plays an important role for the echo signal. Numerical calcu- 
lations of this quantity for an atom-optics Bunimovich stadium are presented. 
Next, by overcoming the inhomogeneous shift, which is the main limiting fac- 
tor for the coherence time in dipole traps [8], echo spectroscopy opens the way 
for using dipole trapped atoms in quantum information and possibly for high 
precision spectroscopy. Finally, our analysis leads to a proposal of a dephasing 
free dipole trap. 



2 Principles of Echo Spectroscopy 

One of the most widely used MW spectroscopy techniques is the so called 
Ramsey method [9]. It works on two level atoms after the following principle: 
The atoms are prepared in one of their two states (denoted |1)). Next, a short 
MW-pulse close to resonance transfers the atoms into a coherent superposition 
state of the two levels (This pulse is usually called a | pulse), which is left 
to evolve for a certain time, without MW-irradiation. Due to the detuning 
between the atomic resonance frequency and the MW-frequency, the time 
evolution will cause an atomic phase to develop in the rotating frame of the 
MW field. A second |-pulse maps this phase onto populations of the two levels, 
which can be detected. If the MW-frequency is scanned, this phase difference 
is scanned, and sinusoidal oscillations in the populations, known as Ramsey 
fringes, are seen. When Ramsey spectroscopy is performed on an ensemble of 
inhomogeneously broadened two level atoms, each atom will acquire a slightly 
different phase in the dark period, causing the macroscopic fringe contrast to 
decay when increasing the time between the |-pulses. Thus the Ramsey fringe 
contrast is a measure of the dehpasing of the system. 

The dephasing can be reversed by stimulating a coherence echo. This is done 
by adding a population inverting pulse (a 7r-pulse) between the two |-pulses. 
If the TT-pulse is added exactly in the middle between the two |-pulses each 
part of the superposition state generated by the first |-pulse will have spent 
an equal amount of time in both levels at the time of the second —pulse, 
and will therefore be exactly in phase with each other. This means that a 
coherence echo have appeared at the time of the second |-pulse even for a 
system that have dephased completely before the 7r-pulse. When measuring 
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the populations the coherence echo is therefore observed by seeing that all the 
atoms have returned to the initial state after the |-7r-|-pulse sequence. If the 
time (ti) between the first | and the vr-pulse is kept constant, and the time 
between the vr and the second |-pulse {T2) is swept, the coherence echo is seen 
as a decrease in the population of |2) for T2 = ti (see Fig. 1). The population 
P2 for T2 = Ti is referred in what follows as the "echo signal", where P2 = 
indicates full revival of coherence and -P2 = | indicates complete dephasing. 

Experimental data shown in this paper is taken using two hyperfine levels 
the ground state of ^^Rb atoms trapped in a dipole trap. These two levels 



551/2, F = 2, mp = Oy, denoted and 5Si/2, F = 3, nip = Oy, denoted |2)) 
are separated by the energy difference Ehf = hjjRF where ojhf = 27r x 3.036 
GHz is the hyperfine splitting. For spectroscopy of trapped atoms we must 
consider the entire Hamiltonian including both the internal and the external 
degrees of freedom of the atom. Since the dipole potential is inversely propor- 
tional to the trap laser detuning 5 [10] there is a slightly different potential 
for atoms in different hyperfine statesQ and the Hamiltonian of our trapped 
two level atom can be written as: 



H = Hi\l){l\+H2\2){2\ 

= + ^1 (x)) |1) (1| + + W + EuF^ |2) (2| , (1) 

where p is the atomic center of mass momentum and Vi [V2] the external po- 
tential for an atom in state |1) [|2)], much smaller than E^f- The atoms are 
initially prepared in their internal ground state Their total wave function 
can be written as \E' = ip, where ip represents the motional (external 

degree of freedom) part of their wave function. If a MW field close to reso- 
nance with ojhf is applied, transitions between the eigenstates of the Hamil- 
tonian corresponding to different internal states, can be driven. The small 
size of optical traps (typically on the order of 100 /xm or less) as compared 
to the MW wavelength (~10 cm), ensures that the momentum of the MW 
photon can be neglected (Lamb-Dicke regime [10]). Even though the inter- 
nal states of the atoms are described as a two- level system Eq.l describes 
a multi-level system in two energy manifolds. When the MW-field is added 
to the Hamiltonian in Eq.l and the time-dependent Schrodinger equation is 
formulated it is seen that the transition matrix elements between eigenstates 
of Eq. 1 are given by C„„' = {n' \ n) x Mi^2 where Mi^2 is the free space 
matrix element for the internal state transition, and {n' \ n) is the overlap 
between the initial motional eigenstate of Hi and an eigenstate of H2. There- 
fore, when a strong and short MW pulse is applied, the motional part of the 



^ For these two hyperfine levels the matrix elements for the dipole interaction 
are identical, hence the relative difference between the two optical potentials is 
iyi — V2) /V2 — ^^^^^ ~ 10^'^ for our experimental parameters. 
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Fig. 1. Echo signal (P2) measured as a function of time between 7r-pulse and second 
7r/2-pulse (T2) for fixed ri. A dip in P2 is seen, showing a "coherence echo". T2 = ri 

initial wave function is simply projected into V2 (x). If Vi (x) = V2 (x) then 
clearly (n' | n) = Snn'- For a small enough potential difference (perturbation) 
6V = V2 (x) — Vi (x) the approximation {n \ n') ~ 5nn' is valid. In the general 
case, a projected eigenstate of Hi will not be an eigenstate of H2, and therefore 
will evolve in the new potential, causing the overlap |(n (t = 0) | n (t = r))| to 
decay (|n {t — r)) = exp (^—i^T^ 

If Ramsey spectroscopy is performed on atoms initially prepared in an eigen- 
state of Vi characterized by the quantum number n, then the probability to 
be in the internal state |2) after the 7r/2-7r/2 pulse sequence can be shown to 
be: 

P^^l(l + (n e-<^-™+^")^ cos (A„t)) , (2) 



where oumw is the MW-frequency, r is the time between the pulses, and A„ is 
a generalized, state-dependent, detuning (it reduces to the detuning when only 
two levels are coupled by the MW-field), defined so (n exp (^—i — ujmw + ^n) 



is real and positive. Scanning ujmw for a fixed r yields the usual Ramsey fringes 



n 



with a contrast given by 



Since 



n 



exp 



^■2 



exp (-i - UJMW + ^n) t) 
this contrast serves as a measure of the 
der a small perturbation in the Hamiltonian [11]. 



stability o: 



Anjrj n) = \{n{t = 0) \n{t 

exp {t (f ) r) exp {-z (f ) r) 
: the quantum evolution un- 



r))\. 



n. 



In our experiments atoms are not initially in a single motional eigenstate, but a 
thermal ensemble of atoms incoherently populates many eigenstates. The total 
population in |2) is then given by an average of P2 over the initial thermal 
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ensemble. A„ depends on the initial state, and due to this spread in detuning 
the fringe contrast of the ensemble-averaged P2 decays rapidly even when 
\{n{t = 0) \ n {t = t))\ ~ 1 for all populated states (note the analogy to the 
inhomogeneous broadening of two level atoms). Hence, the decay of the fringe 
contrast in Ramsey spectroscopy of thermal trapped atoms does not teach us 
anything about the dynamics in the trap {\{n {t = 0) \ n {t = t))\) but simply 
measures the spread in A„. This spread also limits the interrogation time for 
high precision spectroscopy. 

Equation 1 indicates that if a special trap can be constructed so that A„ 
is independent of n, then the fast decay of the Ramsey fringes will be sup- 
pressed. An example of such a dephasing-free trap is an inverted pyramid of 
exponentially decaying evanescent waves, in which gravity provides the verti- 
cal confinement. This trap have the special feature that V2 (x) is exactly equal 
to Vi (x) shifted along the vertical axis. Therefore all eigenenergies are shifted 
the same amount and A„ is indeed a fixed number. This makes such a trap an 
excellent candidate for high precision spectroscopy. Note, that in order to use 
traps for this kind of measurements, the suppression of the relative spread in 
A„, should be accompanied by a reduction in their average value, e.g. by in- 
creasing the trap's detuning. Since the proposed trap is a "dark optical trap", 
where atoms are trapped by repulsive light forces mainly in the dark [12], 
hence minimizing the amount of interaction with the trap's light, then in ad- 
dition to being dephasing-free it has also a low photon scattering rate (which 
will now be the limiting factor for the coherence time) and a small A„. 

Apart from using special geometries a possible solution to the above mentioned 
dephasing is to reverse it using "echo spectroscopy" by adding a MW vr-pulse 
between the two | pulses. P2 for atoms initially populating a single eigenstate 
then becomes: 



P2 = - 
2 

_ 1 
~ 2 



1 - Re 



n 



n 



Re e*-^((^„(t = 0) I ipnit = r)) 



(3) 



where {(fn} is a new basis defined by \(pn {t = 0)) = exp (^—i^T^ \n), and 
\(Pn {t = r)) = exp (^—i^rj \(pn {t = 0)). P2 no longer depends on ujmw and 
Ehf and therefore not on any A„ but only on the dynamics^. Ensemble 
averaging over a thermal ensemble therefore does not remove the fingerprints 
of the dynamics. 

We performed echo spectroscopy on ~ 10^ ^^Rb atoms laser cooled to ~ 20 



Re [exp [i^Tj {ipn (t = 0) | ipn {t = T))j ~ \{ipn (t = 0) | ipn (t = t))\ is ofteii a 
good approximation [7], but the limitations of it is not yet known to the authors. 
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Fig. 2. Echo signal P2 measured for three different trap laser wavelengths A close 
to the Di transition wavelength of 795 nm. Lower curve: A "small perturbation" 
(A=805 nm). A good echo signal (P2 << 1/2) is seen almost independent of time 
between pulses. Middle curve: A "medium perturbation" (A=798.25 nm). The echo 
signal oscillates around a value smaller than 1/2 with partial revivals at r = 1/2tosc 
and r = Tqsc-, where rosc=3.6 ms is the measured trap oscillation frequency in the 
transverse direction. Upper curve: A "strong perturbation" (A=796.25 nm). After 
a short time the echo signal completely disappears (P2 = 1/2), but partly revives 
again at r = 3.3 ms, close to Tosc- 



jjiK and then trapped in a ~ 30 iiK deep Gaussian optical trap with waist size 
of ~ 50 jim. The procedure is detailed in [7]. The echo signal P2 was measured 
as a function of time between pulses for different wavelengths of the trap laser, 
thereby changing the perturbation strength 5V ^ while keeping the trap depth 
constant by adjusting the power of the trap beam. For a small perturbation 
a good echo (P2 << 1/2) is seen independent of the time between pulses (see 
Fig. 2). For larger perturbation damped oscillations to a level smaller than 
1/2 are seen. An even larger perturbation causes a complete decay of the echo 
coherence (P2 = 1/2) followed by partial revivals, which are clear signatures 
of the dynamics. Two revivals (of different amplitude) are seen when the time 
between MW pulses is close to \Tosc and Tosc-, respectively, where Tosc = 3.5 
ms is the measured radial oscillation frequency of the trap (the longitudinal 
oscillation frequency is ~ 1 s, hence the longitudinal motion is essentially 
fixed for the duration of our experiment). The revival at \Tosc corresponds to 
wavcpackct oscillations excited by the symmetric change of the trap strength, 
and the stronger revival at T^sc corresponds also to the vertical sloshing of 
wavepackets induced by the vertical shift of the gravito-optical potential for 
1 2) compared to |1) (in our experiments this vertical shift is of the order of a 
few nm, as compared to a trap cloud size of ~ 50 fim\). The lack of complete 
revival at Tosc is mainly due to anharmonicity of the Gaussian trap, but also 
small amount of decoherence due to spontaneous scattering of photons and 
noise may contribute. 
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Fig. 3. Probability density plot of an eigenstate of a Bunimovich stadium (ko ~ 100). 

The sensitivity of our technique to map the quantum djTiamics of a system is 
revealed here, since we see the dynamics due to a perturbation (kick), that is 
about three orders of magnitude smaller than ksT. 



3 Eigenstate Stability and Atom Optic Billiards 

In Fig. 2 it is seen that for sufficient large detuning the long time asymptotic 
behavior of the echo signal from thermal atoms is smaller than 1/2, and even- 
tually becoming ~ 0. This can be understood from Eq. 3: If 6V is small there 
will be, for each initial motional eigenstate |n) of Hi a large matrix element 
for going into the equivalent motional eigenstate of H2, since Hi and H2 are 
almost identical. If the time between pulses is large enough so the dynamics 
has dephased completely between them, then the contributions from other 
states will average to zero. Hence, the long time ensemble average of P2 is 
simply given by P2 = |[1 — \{n' = n|n)|^], where \n) and \n' = n) are corre- 
sponding eigenstates of Hi and H2 respectively. In [7] a paraxial numerical 
calculation of P2; averaged over the ~ 4 x 10^ thermally-populated states of 
the trap yielded good quantitative agreement with measured long time echo 
signal. 

Because of the importance of the eigenstate stability to the revival of coherence 
(good echo signal) and dynamics in the trap, we now calculate | in' = n\n) \ for 
perturbations, that we assume approximate our inherent perturbation, i.e. the 
difference in dipole potential for the two hyperfine states, for different traps. 
Note, that since a large \ {n' = n\n) \ also implies a large survival probability 
as defined in [11,13,14], it also plays an important role in quantum dynamical 
studies. 

Atom optic billiards are good candidates for trapping atoms for high precision 
spectroscopy, since the interaction with the trap light is small because atoms 
are confined mainly in the dark [12]. They have been used to study dynamics in 
the past, and have experimentally been demonstrated to display the features 
of ideal billiards, [15,16,17,18]. For our calculations we therefore approximate 
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Fig. 4. \{n' = n\n)\ as a function of perturbation strength 1 + Sk/k for eigenstates 
with ko ~ 100. Three different perturbations are shown for the stadium eigenstate 
of Fig. 3, namely dilation, stretehing and the "physical perturbation" (a dilation 
followed by a negative stretch). The rectangle state has kx — ky. For dilation sub- 
stantial reduction of | (n' = n|n) | indeed occurs for 5A; ~ 1, as expected from the 
intuitive discussion in the text. 



Fig. 5. I (n' = n\n) \ for three states of the stadium with ko ~ 100 as a function of 
perturbation strength, (a) Stretching: A strong dependence on the initial state is 
seen, (b) Dilation, for the same three stases as in (a): Only a weak dependence on 
the initial state with similar k is seen. 

the atom optics billiard with an ideal hard-wall billiard, and use the power- 
ful technique of Vergini and Saraceno to numerically calculate highly excited 
eigenstates and eigenenergies [19]. For an unpcrtTirbcd Bunimovich stadium 
[20] (two semicircles of radius r = 1 connected by two straight lines of length 
/ = 2, for units defined such that h = 2m = 1) a plot of the probability den- 
sity of a highly excited state is shown in Fig. 3. The size of its wave vector 
is /co — 100 (V^V — "^oV') hence the wavefunction have on the order of few 
thousand antinodes. 

Systems respond very different to different types of perturbations [21,22,23,13,14]. 
Perturbations of Vi (x) that can be described as a change of coordinate system 
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Fig. 6. (a) Matrix elements (| {n'\n) |) between an initial state and its ten closest 
neighboring states with the same parity, as a function of stretching factor, (b) Eigen 
A;'s as a function of stretching factor, for the same states as (a). It is seen that when 
the state undergo avoided crossing it mixes strongly with the crossed state. Bold 
lines in (a) and (b) corresponds to the same states. 

such as rotation, translation (of importance to the inverted pyramid discussed 
above) and dilation will give exceptional high \ {n' — n\n) \. Intuitively this is 
clear since an eigenstate of Hi under these three perturbations is also an eigen- 
state of H2 rotated, translated or dilated respectively. Therefore, to cause a 
significant drop in | (n' = | the perturbation has to move the billiard wall 
distances on the order of typical distances between nodes in the eigenstate 
wave function, corresponding to ~ 1 in case of dilation (see Fig. 4). On the 
other hand perturbations that changes the symmetry of the billiard can cause 
much faster changes in \ {n' = n\n) \. This is illustrated in Fig. 4 where it is 
seen that dilating a Bunimovich stadium causes \ (n' = n\n) \ to drop much 
slower than changing the lengths of the straight sections (this perturbation 
we call "stretching"). We are interested in the perturbation induced in an 
atom-optics billiard by changing the strength of the optical potential. This is 
best approximated by moving the billiard's wall a fixed distance perpendicu- 
lar to itself. This perturbation, that we name "the physical perturbation", can 
be viewed as a dilation followed by a negative stretch, and it shows features 
between those of dilation and stretch (Fig. 4). 
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Fig. 7. k as a function of dilation factor for eigenstates with kQ ~ 100 in a stadium. 
Straight hues are observed as excepted. 

In Fig 4 I (n' = n|n) I for a rectangle is also shown in order to illustrate that the 
criteria stating that the dilation should move the wall of the billiard a distance 
of the order of the distance between nodes is valid not only for chaotic shapes 
but also for regular ones. For dilation in ID systems this criteria is equivalent to 
the necessary criteria for substantial mixing from perturbation theory, namely 
that the perturbation energy should be comparable to the level separation. 
That this criteria also holds for dilating an integrable shape in 2D (where it 
corresponds to a perturbation energy of ~ 100 times the level spacing in our 
case) is also obvious since it is equivalent to two ID systems. But its fulfilment 
for a chaotic system where variables cannot be separated is less obvious unless 
the above arguments are considered. 



n\n) 



as a function of perturbation strength for different 



Figure 5 shows | {n' 

eigenstates with ko — 100. It is seen that | (n' = n\n) \ is highly state dependent 
in case of stretching, but is not in case of dilation. Following the explanation 
given above, it is clear that states with almost the same ko have almost the 
same \ {n' = n\n) \ for a given perturbation strength in case of dilation. The 
sudden drop in \ {n' = n\n) \ in case of stretching is often associated with 
avoided crossings of energy levels, after which the two states involved still hold 
the major part of the probability. This is clearly seen from Fig. 6, where matrix 
elements to the ten closest states of same parity are shown for stretching^. In 
Fig. 6b as a function of stretching factor is shown for the same ten states. As 
seen, the avoided crossing of the two neighboring states (drawn in bold lines) 
occurs for the same stretch factor as the one causing strong mixing between 
them (Fig. 6)). In Fig. 7 A; as a function of the dilation factor is presented, 
showing straight lines as expected. Matrix elements to close by states are 
for dilation, in strong contrast to stretching shown in Fig. 6a [21]. 



^ Stretching is an even parity perturbation and selection rules therefore insures 
coupling only to states with the same parity 
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4 Conclusions 



We have analyzed echo spectroscopy of trapped ultra-cold atoms. Using this 
technique, we can combine the high precision of MW-spectroscopy with the 
high degree of control over experimental parameters provided by laser cooling 
and trapping techniques, to perform quantum dynamical studies even when 
many states are thermally populated. Based on this analysis we proposed a 
gravito-optical trap in which atoms in superposition states does not suffer 
from dephasing due to inhomogeneous Stark shifts. 

We analyzed numerically the stability of eigenstates for several perturbations, 
and verified that the rapid drop in \ {n' = n\n) \ seen for symmetry breaking 
perturbations is linked to avoided crossings of energy levels. The symmetry 
of the effective perturbations, caused by MW transitions in our experimen- 
tal atom optics billiards was shown to depend on their shape. For example, 
rectangle or a circular billiard the physical perturbation corresponds approxi- 
mately to pure dilation, thus enhancing the applicability of echo spectroscopy 
for these shapes, as compared to the Bunimovich stadium billiard discussed 
above. 

Finally, the suppression of dephasing achieved with echo spectroscopy yields 
a dramatic increase in the coherence time for trapped atoms, that may find 
important applications for precision spectroscopy and quantum information 
processing. 
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